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Abstract 

A quantum field theoretical approach, in which a quantum probe is used to investigate 
the properties generic non-flat FLRW space-times is discussed. The probe is identified with a 
conformally coupled massless scalar field defined on a space-time with horizon and the procedure 
to investigate the local properties is realized by the use of Unruh-DeWitt detector and by the 
evaluation of the regularized quantum fiuctuations. In the case of de Sitter space, the coordinate 
independence of our results is checked, and the Gibbons-Hawking temperature is recovered. A 
possible generalization to the electromagnetic probe is also briefly indicated. 

PACS numbers: 04.70.-s, 04.70.Dy 


1 Introduction 

Hawking discovery of black hole radiation [1] is considered one of the most important predictions of 
quantum field theory in curved space-time. The predicted effect is quite robust, see isiEiiaiaiiiiTiiH]. 

Parikh and Wilczek [9] (see also [EI E]) introduced a further approach, the so-called tunnelling 
method, for investigating Hawking radiation. Later, the Hamilton-Jacobi tunnelling method m 
mods] was introduced. This method is covariant and can be generalized to the dynamical case 

na [13 [m [m Ea EI [22] ■ 

The aim of this paper is to continue the investigation of the local properties of a generic FLRW 
space-time by making use of suitable quantum probe, along the line presented in reference [23] . 
We recall that FLRW space-times may be regarded spherically symmetric (dynamical) space-times, 
and a covariant formalism introduced by Hayward is at disposal [181 da EI]. This approach will 
permit the introduction of Kodama vector and related observers. For the sake of completeness, 
first we review this general formalism. 

* rabochaya@science.unitn.it 
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To begin with, let us recall that any spherically symmetric metric can locally be expressed in 
the form 

ds^ = , i,j^ { 0 , 1 } , ( 1 . 1 ) 

where the two-dimensional metric 


= 'yij{x^)dx^dx^ ( 1 . 2 ) 

is referred to as the normal metric, {x*} are associated coordinates and R{x^) is the areal radius, 
considered as a scalar field in the two-dimensional normal space. Let us introduce the scalar 
quantity 

x{x) = {x)diR{x)djR{x). (1.3) 


The dynamical trapping horizon, if it exists, is located in correspondence of 


X{x] 


= 0 . 


H 


(1.4) 


provided that dixln 7^ 0. Hayward surface gravity associated with this dynamical horizon reads 


kh = 


H 


In any spherical symmetric space-time there exists the Kodama vector field JC, defined by 


IC\x) = 


1 


-7 


■.e^^djR, 


/C® = 0 = /C‘^ 


(1.5) 


( 1 . 6 ) 


Kodama observers are characterized by the condition R = Rq. 

Corning back to the Hamilton-Jacobi tunneling method, we recall the semi classical emission 
rate reads 


r oc |Amplitude|^ oc e ^ « 


(1.7) 


with 9 standing for the imaginary part, the appearance of the imaginary part due to the presence 
of the horizon. The leading term in the WKB approximation of the tunnelling probability reads 


r 


oc e "ir 


( 1 . 8 ) 


in which an energy uj of the particle, and the Hayward surface gravity evaluated at horizon appear. 
Here we recall the operational interpretation |18j . We note that static observers in static BH 
space-times become in the dynamical case Kodama observers whose velocity 


Vk = Jijvkvx = -1 • ( 1 - 9 ) 

Kodama observers are such that R = Rq, namely they have constant areal radius, and the energy 
measured by these Kodama observers at fixed Rq is 


E = -vW = 


K^djl _ ^ 


( 1 . 10 ) 
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As a resuit, the tunneling rate can be rewritten as 


r ~ e 




~ e 


E 

'Tq 


( 1 . 11 ) 


and the local quantity Tq at radial radius Rq is invariant, since it contains the invariant factor y/x, 
and the Hayward surface gravity 


To = 


Th 


Th 


kn 

2tt ■ 


( 1 . 12 ) 


In the static case x = = —goo and recalling Tolman’s theorem: for a gravitational system at 

thermal equilibrium, T -^—500 = constant, it follows that Th = ^ is the intrinsic temperature of 
the BH, namely the Hawking temperature. 

As an important example that illustrates the gauge independence of the formalism is the de 
Sitter space-time. We shall consider three patches, or coordinate systems. The first is the static 
patch, namely 


+ 




(1 - Hl, 


+ r^dS^ 


(1.13) 


Here R = r, the horizon is located at rn = and the surface gravity kh = Hq. 

The second patch is the one which describes a exponentially expanding FLRW flat space-time. 


ds^ = -dt^ + (dr^ + r^dS^) . (1.14) 

Here R = the dynamical horizon is Rh = there is no Killing vector, but the Hayward 

surface gravity is again kh = Hq. Finally, there exists the so called global dS patch, a non-flat 
FLRW space-time 

ds^ = -dt^ -k cosh^(ddot) (—— ox + r‘^dS‘^'\ . (1-15) 

V(l-dd^H) J 

Here R = cosh{Hot)r, and a straightforward calculation gives again Rh = and kh = Hq. 

In the dynamical case, but for slow changes in the geometry, the question is: could possibly the 
quantity Th = ^ be interpreted as a dynamical Hawking temperature? In our opinion a complete 
answer is stili missing, see also [22] for a recent discussion. With regard to this issue, it should 
be important to have a quantum field theoretical confirmation of the tunnelling results (see, for 
example m)- 

As a first step toward this aim, we recall that in the flat FLRW case, a conformally coupled 
massless scalar has been used as a quantum probe in order to investigate such space-times (see 
[MI Ea [27] and references therein). The aim of this paper is to extend the investigation to the 
non flat FLRW space-times, since, at least in the case of De Sitter space-time, there exists the 
important example (|1.15l) of non flat FLRW space-time, namely the global de Sitter patch. 

The paper is organized as follows. In Section 2 we give a brief survey of a scalar quantization 
in a generic FLRW space-time. In Section 3, we present the main formula of Unruh-DeWitt 
response function detector. In Section 4, the computation of the renormalized quantum fluctuation 
is presented and the de Sitter case is discussed in detail. In Section 5, the conclusions are reported, 
and in the Appendix A, an elementary derivation of the Wightman function is given. 
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2 Conformal quantum probe in FLRW space-times 

In the following, we review the quantization of a conformally coupled massless scalar field in a 
generic FLRW space-time. First, it is convenient to introduce the conformal time ry by means 
dri = ^. Thus, we have 

= a^(ry)(-dry^ + dS|), , (2.1) 


where the metric of the spatial section may be written 




1 — 


+ r^dS^ 


( 2 . 2 ) 


with A: = 0, ±1, and /iq is a mass or inverse lenght scale, related to the Scalar Ricci curvature, which 
reads 


n” 

7^ = 6 ( ^ + /c-| 
a'* 


a = 


da 

drj 


A useful and equivalent form for the spatial section is 

dT.l = de + hl{OdS\ 

where = sin“^ h^r and 


KOi = 


sin /lo? 
hn 


MOo = e, M0-1 = 


sinh /loC 
hn 


(2.3) 


(2.4) 


(2.5) 


In the case of a free massless scalar field which is conformally coupled to gravity, the related 
Wightman function W{x,x') =< 4>{x)(f>{x') > can be computed in an exact way. For the sake of 
completeness, we re-derive this well known resuit. 

The quantum field (j) has the usual expansion 


Hx) = f»ix)aa + faix)at 

a 


( 2 . 6 ) 


where the modes functions /«(x) satisfy the conformally invariant equation (TZ being the Ricci 
curvature) 


(° ~ ^ ^ ■ 

Defining the vacuum by OqIO >= 0, the Wightman function turns to be 
W(x, x') = ^ fa{x)f*{x '), 

Q; 

and it satisfies, with x' fixed, 
lF(x,x') = 0 . 


(2.7) 


( 2 . 8 ) 


(2.9) 
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An elementary derivation is presented in Appendix A. Here we give another derivation based on 
conformal invariance. 

First we recall that if we make a conformal transformation 


ds'^ = n{xfdsl, (t) = , 


( 2 . 10 ) 


( 2 . 11 ) 


the related Wightman W{x.x') =< (j){x)4>{x') > transforms as 

We are interested in the non flat case. It is sufficient to consider only the k = 1 case. The k = —1 
may be obtained by the substitution /iq —)• ifiQ. 

Recall that the metric on a non flat FLRW positive spatial curvature space-time is 


= a^{ri) ( —dr]'^ + d^‘^ + sin^ ho^dS^ 


hl 


( 2 . 12 ) 


As a resuit, the above metric is conformally related to a static Einstein space R x S3, with metric 

ds'^ = —drf + d^^ + sin^ ho^dS^ , (2-13) 

% 

On the other hand, it is well known that a static Einstein space-time is conformally related to 
Minkowski space-time, since 


V + L 




(2.14) 


(2.15) 


ds% = 4cos^(/io-^-^—)4cos^(/io-^-^—) [—dt^ -|- dr^ -|- 
with the Minkowski coordinates given by 

t±r = ^ tan(/io^^^^). 
tlQ z 

Due to the homogeneity and isotropy of FLRW space-times, we may take W{x,x') = W{x',x) = 
W{rj — r]',r — r'), same radial separation. Thus, since the Minkowski Wightman function is known, 
making use of equations (|2.1ip and (12.1411 . the Einstein space Wightman function turns out to be 

^ (2,16) 


We{x,x') = ^ 


87 r2 cos(ho(ry - ry')) - cos(/io(^ - ^')) ’ 

As a consequence, again making use of ()2.1ip . the Wightman function related to FLRW spherical 
spatial section is given by 

1 


W{x,x') = 


Uq 


(2.17) 


8Tr‘^a{r])a{r/) cos(/io(ry — r]')) — cos(/io(.^ — ^')) 

Finally, the Wightman function related to FLRW hyperbolic spatial section can be obtained by the 
replacement ho —)• iho, and reads 

1 


W{x,x') = — 


hl 


8TT‘^a{rj)a{ri') cosh(ho (?7 — rj')) — cosh(ho(^ — ^')) ' (2-18) 

As a check, the Wightman function related to FLRW flat spatial section can be obtained taking 
the limit /iq —>• 0. Again for same radial separation 

1 1 


W{x,x') = 


4:7r‘^a{rj)a{rf) —(ry — t]')^ + {r — r')^ 

These results are in agreement with the ones obtained in Appendix A. 


(2.19) 
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3 The Unruh-DeWitt detector in non flat FLRW space-times 


The Unruh-DeWitt detector approach is a well known and used technique for exploring quantum 
field theoretical aspects in curved space-time. For a recent review see [2H]. Here, we review the 
basic formula following Ref. [29l |23] . 

The transition probability per unit proper time of the detector depends on the response function 
per unit proper time which, for radial trajectories, at finite time r may be written as Ar = t—tq > 0 

F{E,T) = ^ReJ^ dse-^^^W{T,T- s), (3.1) 


where tq is the detector proper time at which we turn on the detector, and E is the energy associated 
with the excited detector state (we are considering E > 0). The flat case has been already considered 
in several places (see, m), thus we shall consider the k = 1 FLRW case, namely 


FiE,T) 


Re ds 

IGvr^ Jq a{T)a{T — s) cos /io(Ay;(s)) — cos /io(Ar/(s) — iO) 


(3.2) 


where Ax(s) = xi'^) ~ ~ s)-, and A7/(s) = ??(r) — 7/(r — s). The fO-prescription is necessary in 

order to deal with the second order pole at s = 0. However, we will show in the next Section that 
the leading singularity in the coincidence limit, namely for small s, is 

W{t,t - s) c:t + . (3.3) 

As a resuit, one may try to avoid the awkward limit e —>■ O'*' by omitting the e-terms but subtracting 
the leading pole at s = 0 (see [29] for details), and introducing the quantity 


S^(r, s) = a(T)a(T - s)^ (cos /io(Ax(s)) - cos /io(Ar/(s))) , 


hl 


one can present the detector transition probability per unit time in the form 

1 

i;2(r, s) 


1 r 


ds cos{Es) 


+ ^ ) + ’ 


where the ”tail” or finite time fluctuating term is given by 

cos{E s) 


Jr(E) := - 


2vr2 


At S 2 ( r , s ) 


(3.4) 


(3.5) 


(3.6) 


In the important stationary cases (examples are the static black hole, the FLRW de Sitter 
space), one has S(t, s)^ = H^(s) = S^(—s), and Eq. (13.5p simply becomes 


The first term is independent on r, and all the time dependence is contained only in the fluctuating 
tail. 
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3.1 The de Sitter space in FLRW patches 

In order to check the coordinate independence (gauge-invariance) of our resuit, it is instructive to 
investigate the de Sitter space-time. With regard to the other two FLRW patches, they are the 
flat spatial section, physically relevant for inflation and the positive spatial curvature patch. The 
flat case has already been considered, see for example [23]. For the spatial curved patch, we shall 
make use of the formula derived in the previous subsection. 

Recall that in a generic non- flat FLRW space-time, the Kodama observer is given by 

R{t) = Ro = ^ sin hox, (3.8) 

tlQ 

with constant Rq. For a radial trajectory, the proper time in the non-flat FLRW is 


= a?{ri){drf‘ — . 


(3.9) 


Thus, on the Kodama trajectory 




dt^ 


V a^{t)-RlHl)- 


(3.10) 


In the case of de Sitter, we put /iq = Hq and the flat case is simple and one has an explicit 
expression for rj{T) [23]. In the non-flat case, also for de Sitter, it is not easy to get an explicit 
expression of the conformal time as a function of r. For this reason we consider the Rq = 0 case, 
comoving Kodama observer. Thus dr = dt. Furthermore, since a{t) = coshddo^ = cosh Hqt, one 
has 


T/(r) = — arctan . 
Hq 

We have to compute 


(3.11) 


S^(t, s) = -—2 cosh(ddo'r) cosh Hq{t — s) {cos{HoAr]) — 1) 


Hi 


(3.12) 


Making use of well known trigonometric identities, a direct calculation leads to following results 


HoAr]{T, s) 


—2 arctan 


/ sinh \ 

\^coshddo(7' + f)/ ’ 


(3.13) 


sin^ HoAr]{T, s) = 


sinh^ ^ 


cosh^ Ho{t + ^) + sinh^ ’ 


S^(r, s) = —2 cosh(ddo''‘) cosh Hq{t — s) sin^ {HqAij^t, s)) . 

Hq 

As a consequence, the invariant distance (|3.12|] can be re-written in the form 


S^(r, s) = —-^sinh 


Hi 


Hns 


(3.14) 

(3.15) 

(3.16) 
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Since S^(r, s) = S^(s) = S^(—s), we may use (j3.7p and obtain, for > 0 and making use of 
Residue Theorem 


F{E) = 


E 


1 


r)_ 2'kE 

2^ eW - 1 


(3.17) 


which shows that the Unruh-DeWitt detector in the non flat FLRW de Sitter space detects a quan¬ 
tum System in thermal equilibrium at a temperature Tq = ^, Gibbons-Hawking resuit is recovered 
m- This is an important check of the approach, since it shows the coordinate independence of the 
resuit for the important case of de Sitter space. 

The response function for unit proper time, in the stationary cases we have considered, gives 
information about the equilibrium temperature via the Planckian distribution. We also may argue 
as follow. In the stationary case, in the limit r ^ oo, one has 




ds e 


—iEs 


1 


ms) 


Note also that in this case one has 
F{E) 




1 ^ 1 


E 


E 


exp 


To 


- 1 


(3.18) 


(3.19) 


Viceversa, if the above relations hold then E{E) is the Plank distribution. Thus we may define the 
local equilibrium temperature by means of 


Tn = 


E 


\nE{-E) -lnE{E) 


or 


Tn = 


E 


In (l + — 


(3.20) 


(3.21) 


2n{F)J 

which shows in which sense a Unruh-DeWitt detector is a quantum thermometer. 


3.2 d-dimensional generalization 

In the de Sitter case, it is possible to generalize the computation of the response function per unit 
time to the massive non conformally coupled d-dimensional scalar field m- 
In this case one may directly obtain 


FdAE) = 


, n ZL® 




d — 1 v 




E , ,d — 1 V 


. E , 

-h l - 

4 2 2Ho' 


(3.22) 


where 


' (d — 1)2 m2 


u = 


(3.23) 
















with m? mass of the scalar field and ^ the coupling constant with the Ricci curvature, being the 
conformal coupling ^ consequence, one has 

Fd,u{-E) = PdAE). (3.24) 

Furthermore, the massless conformally coupled case in d dimension corresponds to = ^. In this 
case, for d = 4 one recovers 


E^l{E) = — 


E 


1 


r)_ 2ttE 

2vr 1 


Furthermore, making use of the identity 
|2|t-/i m 2 27r2|z|2 


riz)mi-zr = 


cosh(27rImz) — cos(27rRez) 


one also has for d = 3 and ^ 



1 / E^\ 1 

16 e ^+1 


(3.25) 


(3.26) 


(3.27) 


in which the well known phenomenon of the inversion of the statistic for odd dimensional spaces is 
manifest. 

The other interesting case is the minimally coupled massless scalar field, for which u = In 
particular for d = 4, one has 





(3.28) 


In this physical important case, it mimics the graviton, the well known infrared problem associated 
with it appears in the bad behavior for small E. 


4 Quantum fluctuations 

Another proposal to detect local temperature associated with stationary space-time admitting an 
event horizon has been put forward by Buchholz and collaborators [32] (see also |33ll3l|). The idea 
may be substantiated by the following argument. 

Let us start with a free massless quantum scalar field (p{x) in thermal equilibrium at temperature 
T in flat space-time. It is well known that finite temperature field theory effects of this kind may 
be investigated by considering the scalar field defined in the Euclidean manifold Si x where 
one has introduced the imaginary time r = —R, compactified in the circle S*!, with period /3 = y 
(see for example [55]i. 

Let us consider the local quantity < ())(x)^ >, the quantum fluctuation. Formally this is a 
divergent quantity, since one is dealing with product of a valued operator distribution in the same 
point X, and regularization and renormalization are necessary. A simple and powerful way to 
deal with a regularized quantity is to make use of zeta-function regularization (see for example 
[MI [371135], and references therein). Within zeta-function regularization, one has 

<(j){xf >=C{l\Lp){x), (4.1) 
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where (^{z\Lp){x) is the analytic continuatiori of the local zeta-function associated with the operator 

Lp 

Lp = -dl - V2 , (4.2) 

defined on Si x E?. The local zeta-function is defined with Rez sufhciently large by means 
1 

C,{z\Lp){x) =-— j dtE~^Kt{x,x), (4.3) 

where the heat-kernel on the diagonal is given by 

1 


Kt{x,x)=<x\e 


/3(47rt)3/2 


e ^ 


(4.4) 


In (14.ip the analytic continuation of the local C.i^\Lp)ix) appears and it is assumed that this 
analytical continuation is regular at 2 = 1, which, as we shall see, is our case. If the analytic 
continuation has a simple pole in z = 1, the prescription has to be modified (see [38l|33]). 

A Standard computation, which makes use of the Jacobi-Poisson formula leads to 


Kt{x,x) 


1 

(47rt)2 



(4.5) 


Let us plug this expression in (|4.3I) . The term n = 0 leads to a formally divergent integral 
but this is zero in the sense of Gelfand analytic continuation, and it can be neglected. 
Thus, a direct computation gives the analytic continuation of the local zeta-function 

(t) *' 

where Cr{z) is the Riemann zeta-function. It is easy to see that the analytic continuation of the 
local zeta-function is regular ai z = 1, and froin (14.ip . recalling that Cj?( 2) = one has 


< 4>{x)‘^ >= 


1 

12/32 


T2 
12 ■ 


(4.7) 


Thus, the zeta-function renormalized vacuum expectation value of the observable (fP', the fluctua- 
tion, gives the temperature of the quantum field in thermal equilibrium, namely one is dealing with 
a quantum thermometer. 

Motivated by this argument, let us consider again a conformal coupled scalar field in a FLRW 
non flat space-time. We have seen that the off-diagonal Wigthman function is 


W{x,x') =< 4’{x)4>{x') >= 


1 1 
47r2 ^ 2 ( 3 ;^ x') ’ 


(4.8) 


where 


S^(r,r - s) = a{T)a(T - s)^ (-cos Fo(AC(s)) + cos , 

-^0 

with a(r) being the conformal factor. In the limit s —?• 0, formally one has 


(4.9) 


< (/>(x)^ >= W{t,t) , 


(4.10) 
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but W{t,t) is ili defined, and one has to regularize and then renormalize this object. In our case, 
we may make use of the simple point splitting regularization [3], namely we consider W{t,t — s) 
and evaluate the limit s —>■ 0. 

To implement this limit procedure, one has to make use of several identities. For radial time-like 
separation, the starting point is 

«^(^) -X^) = 1- (4-11) 


Taking first and second derivativos with respect the proper time, one has 

^ + = 0, (4.12) 

and 

aHnV - a = f) - - + 3(-)2 . (4.13) 

a a 

Making use of these identities, a long but straightforward calculation leads to 

= -s2 ^ _ (4 

where 

h'^ 

B = + 2Hi + (1 - 2f) . (4.15) 

In this expression A? is the square of the 4-acceleration along the time-like trajectory, given by 
A^ = ^^{t + H{t^-l)) . (4.16) 


In the above expression, the last term is the new one with respect the flat case. Thus, the point 
splitting gives 

With regard to the renormalization, we simply subtract the first divergent term for s —>■ 0. The 
physical meaning of this subtraction has been discussed in detail in reference [23], and it amounts 
to subtract the contribution related to an inertial trajectory in Minkowski space-time. Thus, the 
renormalized quantum fluctuation reads 

< >R= + A‘^ + 2Hi + ^{l- 2t2)^ . (4.18) 


This resuit is the generalization of the one obtained in a flat FLRW space-time in [23] and within 
Unruh-de Witt detector in [26] . 

As a first important check, let us consider again the de Sitter space-time in the global patch. 
In this stationary case, the fluctuation acts as a quantum thermometer. Recall that in this case, 
one has ho = Hq, 


a{t) = cosh Hot , 


H{t) = Ho 


sin hHot 
cosh Hot ’ 


H = i 


_JA_ 

(coshFfo^)^ 


(4.19) 
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As a consequence 

<4>^>R=^{hI + A^) . (4.20) 

This is in agreement with the resuit obtained for the de Sitter space-time in reference [39j. In fact, 
the acceleration can be computed, since, for a Kodama observer R = Rq, one has 


f-l = 


1 - 


Thus 


:■ -^0 rr(,.dH 

Taking equation (j4.16l) into account, one gets 


(4.21) 


(4.22) 


A2 


1 - RlHl ’ 


(4.23) 


which coincides with the Kodama de Sitter acceleration evaluated in the flat patch |23] . Furthermore 
we also have 


< >R= 


1 

487r2 1 - R^H^ ■ 


(4.24) 


For a comoving Kodama observer Rq = 0, and one gets again Gibbons-Hawking temperature 
associated with de Sitter space-time. 


5 Conclusion 


In this paper, with the aim to better understand the temperature-versus-surface gravity paradigm, 
the asymptotic results obtained by tunnelling semi classical method have been tested with quan¬ 
tum field theory techniques like the Unruh-DeWitt detector and the evaluation of the quantum 
fluctuation associated with a massless conformally coupled scalar field. More precisely the results 
obtained in reference [23] have been extended to a generic non spatial flat FLRW space-time. The 
de Sitter space-time in the global non flat FLRW patch has been used as important example, and 
the Gibbons-hawking temperature for the de Sitter space has been re-derived with our general 
fluctuation formula as well as the Unruh-de Witt detector technique. 

With regard to possible generalization, at least in the flat case, our approach may also be 
extended to the Maxwell field as soon as one makes use of the resuit obtained in reference [40|. In 
fact there, quantizing the Maxwell field in a flat FLRW space-time and in the so called W gauge, 
a conformal lifting of the Lorenz gauge in Minkowski space-time, the Maxwell Wightman function 
has been obtained in the form 


^ Y dyX J Q/{X j J 


(5.1) 


where W{x,x') is the Whightman associated with massless conformally coupled scalar field in 
flat FLRW space-time. This strongly suggests that our conformally coupled scalar probe may 
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mimic quite well the quantum Maxwell field. Of course, working with the Maxwell field, the gauge 
invariant has to be implemented and the relevant quantity is the Wightman function associated 
to (for example) the magnetic field. This is a very interesting issue with important cosmological 
applications, see for example the recent discussion appeared in [HI 02], and we hope to consider 
this issue in a future work. 


Acknowledgement 

We thank G. Cognola, M. Rinaldi, L. Vanzo for several discussions. 


6 Appendix A 

In this Appendix an elementary derivation of Wightman function for a massless conformally coupled 
scalar field is presented. Due to the homogeneity and isotropy of FLRW space-times, we may take 
W (x, x') = W {x', x) = W{ri — 7]' ,r — r'). 

It is convenient to introduce the auxiliary quantity 


Y{x, x') = a{r])a{r]')W(x, x'), 

and, from the equation of motion in FLRW space-time with the conformal time one has 
cPY 


drj'^ 


- kH^Y + AhY = 0, 


( 6 . 1 ) 


( 6 . 2 ) 


where is the Laplace operator associated with the metric dSg. 

We may take x' = 0. Let us start with the flat case k = 0. One has 

q2y QY _ 

drf dr'^ clr 

The solution is 


Y = 


1 


— 7 ]^ _|_ 


(6.3) 


(6.4) 


As a resuit, making use of the homogeneity and isotropy, and dealing with the distribution nature 
of VF, 


W{x,x') = 


1 


1 


47 r^a(r7)a(r/') (r — r')^ — (r/ — r/' — ie)^ 


Here we leave understood the limit as e —>• 0'*'. 

We may rewrite it in covariant form, according to Takagi 
Schlicht’s proposal to our conformally flat case, namely 


and Schlicht 


W{x,x') = 


1 


1 


47 r2a(r/)a(r/') [(x — x') — ie{x + x')Y 
where an over dot stands for derivative with respect to proper time (see also [23]). 


(6.5) 

We adapt 

( 6 . 6 ) 
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Corning back to the non-flat case, it is sufficient to consider the positive curvature case /iq > 0. 
The negative one may be obtained with the replacement —>• ih^. 

For the sake of simplicity, here we may take /iq = 1, then dimensional analysis will give the complete 
expression. Let us start with 


ds^ = a^{r]) {—drf + dx^ + sin^ xdS"^) 
The Ricci scalar reads 


7 ^ = 6 


a" + A 


da 


a = 


a*" dri 


and the equation for Y is 

d‘^Y d‘^Y cosxdY 

— -^ + -^—2 + = 0 • 

^X siny ox 

The solution of this partial differential equation is 
1 


Y = 


cos r] — cos X 

Making use of dimensional analysis one arrives at 

1 


W{x,x') = 


hl 


87 r‘^a{r])a{r/) cos /io(x “ X') ~ cos — 1 ]' — ie) 


(6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 


In the case of negative spatial curvature one has 


W{x,x') 


1 _ ^ _ 

8 TT‘^a{r])a{r]') cosh/io(x — x') ~ coshho{r] — r]' — ie) 


( 6 . 12 ) 
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